arXiv:astro-ph/9701219vl 28 Jan 1997 


To appear in the Astrophysical Journal 


THE BEAMING PATTERN OF DOPPLER BOOSTED 
THERMAL ANNIHILATION RADIATION: APPLICATION TO 

MEV BLAZARS 

Jeffrey G. Skibo and Charles D. Dernier 

E. O. Hulburt Center for Space Research, Code 7653, Naval Research Laboratory, 

Washington, DC 20375-5352 

and 

Reinhard Schlickeiser 

Max-Planck-Institut fiir Radioastronomic, Postfach 20 24, 53010 Bonn, Germany 

ABSTRACT 

The beaming pattern of thermal annihilation radiation is broader than the 
beaming pattern produced by isotropic nonthermal electrons and positrons in 
the jets of radio-emitting active galactic nuclei which Compton scatter photons 
from an external isotropic radiation field. Thus blueshifted thermal annihilation 
radiation can provide the dominant contribution to the high-energy radiation 
spectrum at observing angles 9 ^ 1/T, where T is the bulk Lorentz factor of the 
outflowing plasma. This effect may account for the spectral features of MeV 
blazars discovered with the Compton Telescope on the Compton Gamma Ray 
Observatory. Coordinated gamma-ray observations of annihilation line radiation 
to infer Doppler factors and VLBI radio observations to measure transverse 
angular speeds of outflowing plasma blobs can be used to determine the Hubble 
constant. 


Subject headings: active galactic nuclei — cosmology: distance scale — galaxies: 
radio — gamma rays: galaxies — gamma rays: theory 


1. Introduction 


Observations made with the Compton Telescope (Comptel) on the Compton Gamma- 
Ray Observatory ( CGRO ) reveal the existence of a class of gamma-ray blazars displaying 
spectral energy distributions peaking in a rather narrow energy band centered at a few 
MeV (Bloemen et al. 1995; Blom et al. 1995). By combining data from 19 September - 3 
October 1991 and 27 December 1991 - 10 January 1992, Bloemen et ah (1995) discovered 
the first “MeV blazar” GRO J0516-609, which shows a 3-4 o peak in its vF v spectrum in 
the 1-10 MeV energy interval. Only upper limits are measured with Comptel in the 0.75-1 
MeV and 10-30 MeV range, although contemporaneous observations with the Energetic 
Gamma Ray Experiment Telescope (EGRET) on CGRO indicate the presence of a weak 
3 cr excess at photon energies E > 100 MeV. The flat spectrum quasars PKS 0506-612 
and PKS 0522-611 are both within the 95% Comptel contours of GRO J0516-609, but the 
former source is the preferred identification due to the more precise EGRET localization, 
which places only PKS 0506-612 within its 95% confidence contour. This source is reported 
at the 3.9a significance level in the Second EGRET catalog (Thomson et al. 1995), but is 
too weak to be listed as a catalogued EGRET source. 

The blazar PKS 0208-512 was detected in analysis of Comptel data by Blom et al. 
(1995). The strongest signal from PKS 0208-512 was obtained by combining data from 8-13 
May 1993 and 3-14 June 1993, yielding a flux in the 1-3 MeV band of 4.1(±0.7) x 10 -4 
photon cm -2 s ^ 1 and upper limits at lower and higher energies in the Comptel energy range. 
Contemporary observations with EGRET yield a strong >100 MeV detection of the source 
2EG J0210-5051, with PKS 0208-512 located within the 95% confidence contour of this 
source. The peak in the uF u spectrum of PKS 0208-512 occurs at MeV energies, supporting 
its classification as an MeV blazar. Williams et al. (1995) also report the detection of 
the gamma-ray source GRO J1753+57 with a spectrum similar to that of the other two 
MeV blazars. Lack of association of this object with a flat radio spectrum quasar makes 
this object’s classification uncertain, but one should keep in mind that blazars are highly 
variable at all wavelengths, so that contemporaneous multiwavelength campaigns could be 
necessary for proper source identification. 

The peaked emission from MeV blazars has been interpreted as Doppler boosted e + e - 
annihilation radiation produced in collimated jets of plasma ejected from a supermassive 
black hole (Roland & Hermsen 1995; Henri, Pelletier, & Roland 1993). The possible 
importance of e + e _ processes and 7-7 opacity effects has been considered by Blandford 
& Levinson (1994), Bottcher, Mause, & Schlickeiser (1997), and Marcowith, Henri, & 
Pelletier (1995). No attempt has been made, however, to account for the significant spectral 
differences between MeV blazars and the > 50 7 -ray blazars detected by EGRET (e.g., von 
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Montigny et al. 1995), which show strong high-energy radiation extending to GeV and, in 
a few cases (e.g., Macomb et al. 1996; Schubnell et al. 1996), TeV energies. 

In this paper we show that the beaming pattern of thermal Doppler-boosted e + e - 
annihilation radiation is much broader than the beaming pattern produced by nonthermal 
jet electrons which Compton scatter photons produced outside the jet, whose beaming 
pattern was recently derived (Dernier 1995; Dernier, Sturner, & Schlickeiser 1997). We call 
this latter process external Compton scattering (ECS) and distinguish it from synchrotron 
self-Compton (SSC) scattering, which has also been proposed as the radiation process 
yielding blazar gamma radiation (see Marscher & Travis 1996 and references therein). 

If the high-energy emission in blazars is primarily produced by ECS and, moreover, the 
annihilation luminosity from thermal e + e” pairs is significant, then a class of objects with 
properties similar to the MeV blazars is a natural consequence of orientation effects when 
viewing at moderate angles with respect to the axis of the outflowing jet. In addition, we 
show in Appendix A that if the peaks of MeV blazars are indeed due to Doppler-boosted 
thermal annihilation radiation, then gamma-ray observations to measure Doppler factors 
T> coordinated with VLBI radio observations to measure transverse angular speeds fi of 
outflowing radio blobs in two or more episodes is sufficient to determine the Hubble constant 
Ho- 

2. Spectra and Angular Distribution of Jet and Accretion-Disk Radiation 

We calculate the received spectral flux and radiation pattern from a system consisting 
of an accretion disk and a relativistic collimated outflow of thermal and nonthermal e + e” 
pairs. We do not attempt to derive the particle distributions self-consistently from first 
principles. Evolution of nonthermal particles due to radiative losses in outflowing plasma 
blobs has been treated in detail by Dernier & Schlickeiser (1993) and Dernier et al. (1997), 
and the formation of a thermalized pair plasma in jets has been considered by Bottcher & 
Schlickeiser (1996). The emphasis here is rather on the angle-dependence of the received 
fluxes from the following three spectral components: (1) annihilation radiation from hot 
thermal pairs in the relativistic outflow; (2) inverse Compton emission from nonthermal 
pairs in the outflow; and (3) a quasi-isotropic thermal Comptonization spectrum from 
hot optically thin accretion plasma located near the central supermassive black hole. The 
central assumption in calculating the jet spectra, as discussed by Dernier & Schlickeiser 
(1993), is that the particle distribution function is isotropic in the comoving blob frame. 

We derive the received spectrum of thermal annihilation radiation. The observed flux 
density S(c, /i, (ft) (erg cm -2 s -1 e _1 ) from a blob moving with bulk velocity /3c is related 
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to the emissivity in the comoving frame j*(e*, n*, <fi*) (erg cm 3 s 1 e 1 sr x ) through the 
relation 


S(c , fi, (/)) 


V\l + z)V* 

j 2 ■) 


e(l + z) 

V 


~P 

-/V 


( 1 ) 


(e.g. Begelman, Blandford, & Rees 1984), where e = hu/m e c 2 is the dimensionless 
photon energy in the stationary frame and V* is the volume of the blob in the comoving 
frame. (Hereafter starred quantities refer to the comoving frame, unstarred quantities 
to the stationary frame.) The direction from the source to the observer is (/i, 0), with 
ju representing the cosine of the angle between the jet axis and the line of sight, and (j) 
representing the azimuthal angle. The Doppler factor T> = [T(l — (3/1 )] _1 , where the bulk 
Lorentz factor T = (1 — /3 2 ) -1 / 2 , the speed of the blob is /3c, and the luminosity distance 
<7 l = 2c[l + z — (1 + z) 1 ! 2 ]/Hq for a 8ttGp/3Hq = 1, A = 0 cosmology (Weinberg 1972). The 
Hubble constant Hq = 50 h km s ^ 1 Mpc -1 . 

We assume that there is a population of thermal pairs with dimensionless temperature 
© a = kT/m e c 2 in the outflow. The annihilation emissivity is isotropic in the comoving 
frame and is given by the expression 


j:(e*,/T, </>*;©„ 
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(Svensson 1983; Dermer 1984), where L* is the luminosity of annihilation photons in the 
comoving frame and R '2 is the modified Bessel function of the second kind. The integration 
variable y r = U + ■ is the invariant relative Lorentz factor formed from the e + and e - 
four velocities. The cross section for annihilation, expressed in terms of y r , is given by 


C r a( 7r) 


(e.g. Jauch & Rohrlich 1976), 
of the electron. 


7n~e jV 7 2 + 47 r + l \ 

7 r + llV 7 2 -l / 

ln[7r + (ifr - 1) V2 ] - (7 J r _ + ^i/ 2 } ( 3 ) 

where r e = e 2 /m e c 2 = 2.82 x 10~ 13 cm is the classical radius 


Combining equations (1), (2), and (3) we obtain the annihilation flux density 
S a (c, The beaming pattern for thermal annihilation radiation is not obvious from 

the resulting formula. To derive the beaming pattern we approximate the annihilation 
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emissivity as a delta function centered at unity, which is a reasonable approximation when 
0 a €l. Thus 

!)■ ( 4 ) 


Inserting this into equation (1) and integrating over e, we find that the annihilation flux 
density integrated over photon energy has a beaming pattern S a oc T> 4 . 

Electron-positron thermal bremsstrahlung radiation will always accompany the 
annihilation radiation. At temperatures 0 a 3> 1, the luminosity of e + -e~ bremsstrahlung 
can dominate e + -e _ annihilation radiation (see Skibo et al. 1995). For the temperatures 
we use in this paper, however, the bremsstrahlung flux is less than the fluxes of the other 
components, and is not included in the fits. 

In addition to the annihilation radiation, we assume that there is an isotropic 
population of nonthermal pairs in the comoving frame with an energy spectrum given by a 
broken power law, so that 

-tat™*) 11 TV" ■ (5) 

4vr ( (7/7 6 ) v 7fe < 7 < 72 

where p 2 = pi + 1 and pi ^ 2. Such a time-averaged particle spectrum might arise if 
shocks are present and radiative cooling is efficient for pairs in the high energy protion of 
the spectrum (Dernier & Schlickeiser 1993; Sikora, Begelman, & Rees 1994). We consider 
a situation where the particles scatter external ambient soft photons which are isotropic in 
the frame through which the radiating plasma moves (Sikora et al. 1994). If the photon 
energy is 6 q, then all scattering will take place in the Thomson regime if Ty^g <C 1. If the 
scattering takes place in the Thomson regime, then the approximation of Reynolds (1982) 
can be used to derive (see Dernier 1995; Dernier et al. 1997) the observed scattered flux 
density, given by 
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for i = 1, 2, where a* = ( p.i — l)/2. When i — 1, 


and when i = 2, 
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The term L* c in equation (6) represents the luminosity of Compton-scattered radiation 
measured in the comoving plasma frame. The factor A in the normalization of expression 
(6) depends only on kinematic quanties through the expression 


A = 


167rT 2 eo 


pi/ 2—2ai 2— 2 qi ' 

% (7 b 


— 7i 
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The beaming pattern in equation ( 6 ) goes as V l+2a {l + p) 1+a ~ T )4+2 " f or sca tt e red 
radiation. If p\ ~ 2, then the spectral index a ~ 1/2 at energies below the break energy, 
and S c (e) oc D 5 ( 1 + yu) 3//2 . Thus the scattered radiation is more highly beamed than 
the annihilation radiation, where S a oc T> A . We note that if the emission is made by a 
radiation process which produces a photon source which is isotropic in the comoving frame, 
then S'(e) oc V 3+a . Such beaming patterns result from isotropically distributed electrons 
radiating synchrotron and synchrotron self-Compton emission in a spherical plasma blob 
(e.g., Marscher 1987). This gives a beaming pattern comparable to the annihilation 
radiation, and when such emission components dominate, we expect the ratios of the flux 
densities of synchrotron or synchrotron self-Compton and annihilation radiation to vary 
only weakly with observing angle. 

The third spectral component we consider is the accretion disk radiation. If optically 
thin, this component is isotropic in the stationary frame and would only display an angular 
dependence if obscuring clouds or a dusty torus shadow the radiation. Such emission could 
arise from Comptonization of soft accretion disk photons by a hot disk corona (e.g., Haardt 
& Maraschi 1993). For simplicity, we model this component by an exponentially truncated 
power law. The observed flux density is then given by 

L d ( 1 + zy-° d e~° d exp [-(1 + z)e/Q d ] 
d[ ) dvr dl dec-**-exp (-£/©„)’ 1 J 

where the disk luminosity is denoted by L d . 

In Fig. 1 we show spectra for various viewing angles, choosing T = 5. The pair 
temperature in the outflow is fixed at © a = 0.5 and the parameters of the nonthermal 
pair distribution are 71 = 1, 7 *, = 10 3 , 72 = 10 6 , pi = 2.1 and P 2 = 3.1. We set a d = 0.5, 
Q d = 0.5, and e c = 10 ~ 3 for the disk parameters and set the redshift z — 1. The disk 
luminosity L d is assigned the value 10 4 ' ergs s _1 . We set L* c = L d /T. The reason for this 
choice is that the jet luminosity should normally be less than the disk luminosity, since 
the accretion power is supplying the jet luminosity, and the factor accounts for the 
additional power required by the bulk motion of the jet. For the annihilation luminosity, we 
take L* = 0AL* c . 
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As can be seen in Fig. 1, the disk radiation is independent of observing angle and 
the annihilation radiation is less beamed than the strongly anisotropic Compton-scattered 
radiation. For this simple model, the disk produces X-ray emission only, the Compton 
scattered radiation always dominates at gamma-ray energies above 10 MeV and the 
annihilation radiation appears as a enhancement in the MeV range which, for certain 
observing angles, can dominate the flux at MeV energies. Hence for small observing angles, 
the source appears as a gamma-ray emitting blazar with the luminosity dominated by 
high-energy gamma-rays and the source displays properties similar to a large fraction of 
the > 50 gamma-ray blazars detected with EGRET (see, e.g., von Montigny et al. 1995). 
At angles slightly off axis the high-energy radiation flux decreases and the bulk of the 
high-energy gamma-ray luminosity comes out in the MeV range. The annihilation feature 
can dominate at moderate observing angles. At larger angles, the disk radiation dominates 
the observed flux, although scattered jet radiation could make a substantial contribution 
to the flux at MeV energies (Skibo, Dernier, & Kinzer 1994). Note also that radio loud 
AGN could exhibit time-variable hardenings and features at MeV energies due to the 
electron-positron component, as suggested in the 6 = 30° panel in Fig. 1. 

We use this model to fit the Comptel data for the blazars PKS 0208-512 at redshift 
z = 1.003 and PKS 0506-612 at z = 1.093 (Bloemen et al. 1995; Blorn et al. 1995). In Figs. 
2 and 3 we show the contemporaneous gamma-ray observations of these sources made with 
the Comptel and EGRET instruments on CGRO. The X-ray points are noncontemporaneous 
observations made with ROSAT (Brinkman et al. 1994). In these fits, we fix the relative 
luminosities of the various spectral components to satisfy L* a = Lq = Y~ 1 L d . For the case 
of PKS 0208-512, the bulk Lorentz factor Y = 3, the inclination angle 6 = 18° and the 
luminosity L d = 10 47 erg s _1 . For PKS 0506-612, we let Y = 11, 9 = 9° and L d = 3 x 10 46 
erg s -1 . All other parameters are the same as used for Fig. 1. The parameters used in the 
fits are not unique, but are typical for gamma-ray blazars. Thus if jets are composed of 
electrons and positrons which thermally annihilate, then the existence of a class of MeV 
blazars could be a straightforward consequence of the different beaming patterns of the 
radiations. 


3. Discussion 

In this paper, we have considered implications of the suggestion that the features in 
MeV blazars arise from thermal annihilation radiation in a relativistic jet, as suggested by 
previous authors (see Section 1). The existence of a class of MeV blazars can be understood 
by orientation effects if the beaming patterns of the nonthermal radiation is different from 



that of the annihilation radiation. If the gamma rays are produced by nonthermal jet 
electrons scattering photons produced external to the jet, then the annihilation radiation 
has a broader beaming pattern. This could account for the existence of the class of MeV 
blazars, and would be consistent with unification scenarios for radio-emitting AGNs (see, 
e.g., Urry & Padovani 1995 for a recent review). Moreover, if this interpretation is correct, 
then we predict that radio galaxies will exhibit time variable hardenings and features at 
MeV energies which would correlate with increasing core dominance. 

Gamma-ray observations of line features and contemporaneous VLBI measurements 
of transverse angular speeds from relativistic plasma outflows provide a method to 
determine the Hubble constant. The radio observations furnish the apparent speeds, and 
the gamma-ray observations give the Doppler and cosmologically shifted energy of the 
annihilation line. As noted earlier by Roland & Hermsen (1995), measurements of V 
and the apparent superluminal speed can be used to determine the angle of the jet with 
respect to the observer and the bulk Lorentz factor T of the outflowing plasma. As shown 
in Appendix A, multiple contemporaneous observations of the shifted line energies and 
transverse angular speeds provide a method to determine the Hubble constant. This is an 
extension of the work of Schlickeiser & Dernier (1995), who proposed a model-dependent 
method for measuring V from the broadband spectral energy distribution of blazars. If 
the MeV features can be conclusively identified with annihilation line signatures, then the 
procedure proposed here will provide a more definitive method for determining the Hubble 
constant from coordinated radio and gamma ray observations. 

Arbitrary spectral lines emitted from the outflowing plasma jets can also be utilized by 
this method. The identification of spectral line features with a specific outflowing plasma 
blob is indicated at gamma-ray energies by time variability. Observations (Wehrle 1996) 
show that the emergence of a radio emitting feature in high-resolution VLBI maps is 
preceded by a gamma-ray flare. A better method would be through milliarcsecond imaging 
of spectral lines, which is presently only feasible at radio and optical wavelengths; high 
resolution optical imaging requires, however, bright sources of apparent magnitudes m <; 9. 

If it is established that the peaks of MeV blazars are a consequence of electron-positron 
plasma jets, then such a result would have important implications for processes which power 
jets in radio-loud AGNs, requiring energization of particles in the compact environment 
near the supermassive black hole. Additional measurements with Comptcl on CGRO, 
and observations with more sensitive soft gamma-ray telescopes, such as the INTEGRAL 
telescope scheduled for launch early in the next millenium, could resolve the question of the 
composition of jets in AGNs and provide a test of this model. 
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A. A Method for Determining the Hubble Constant 


Here we show how gamma-ray observations may be used in conjunction with VLB1 
rdio observations to measure the Hubble constant. Consider a source at redshift z which 
ejects two e + e~ pair plasma blobs with velocities f3\C and P 2 c at different times in the same 
direction # measured with respect to the line of sight. Let <5i and S 2 be the blueshifts of the 
annihilation line emitted by these blobs measured with y-ray observations. We have 

= [Pi(l — Pi cos#)(l + z)\~\ (Al) 

£2 = [r 2 (l — P 2 cos#)(l + z)\ 1 , (A2) 

where T* = (1 — Pf)~^ 1 i — 1,2. Eqs. (Al) and (A2) assume that the peak of the 
annihilation radiation occurs at m e c 2 , whereas the peak of thermal annihilation radiation is 
blueshifted (Ramaty & Meszaros 1981). The blueshifting of the peak of the radiation can 
be accounted for by determining the temperature through the width of the emission feature. 
If this effect is properly taken into account, then 8\ and <5 2 can be precisely determined. 


Let pi and p 2 be the angular speeds of the two blobs determined from radio 
observations. The apparent speeds of the blobs in the plane of the sky are 


c 

d-M l-l 2 
C 


Pi sin 6 
1 — Pi cos 6 
P 2 sin 6 
1 — p 2 cos 6 


= (1 + z)PiTi5i sin#, 
= (1 + z)P 2 T 2 6 2 sin#. 


(A3) 

(A4) 


where cIm is the proper motion distance. Setting the cosmological deceleration parameter 
g 0 = 1/2 the proper motion distance is a function of the cosmological redshift z and the 
Hubble constant Hq (Weinberg 1972) through the expression 


du — 


2 c r 
H n 


(! + *)■ 
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These five equations express the five unknown quanties g?m, #, Pi, P 2 and H 0 in terms 
of the measured quantities pi, p 2 , 5i, S 2 and z. Solving equations (Al) and (A2) for cos# 
and equations (A3) and (A4) for sin #, and equating the ratio tan # to eliminate #, we obtain 
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Dividing equation (A3) by (A4) leads to 


r?=l+ (St) 

Substituting (16) into (17) we arrive at the following expression for Y { : 


Ti = 


---( 1 — — 

Ai (1 + z) y ii 2 , 


-i 


9 1^2 
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fr-g) 
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From equation (Al) we then determine 0, given by 

e = cos -1 !/?] -1 - [/Ar^i + z)] -1 }. 

Using equation (A3) for the proper motion, we obtain 

c (r? - if / 2 


dhi — 


Hi 


-(1 + z)5i sin 9, 


and finally from equation (A5) we find the Hubble constant, given by 

2 c r 


Hn = 
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1 - (1 + z)- 1/2 
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(All) 
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Fig. 1.— Spectra for various viewing angles and bulk Lorentz factor V = 5. The dashed 
curves represent the observed flux produced by inverse Compton scattering of isotropic soft 
photons by nonthermal pairs in the relativistic outflowing plasma jet, the dotted curves are 
the observed fluxes of thermal pair annihilation radiation produced in the jet, and the dot- 
dashed curves are the optically thin isotropic emission produced by an accretion disk corona. 
The solid curves give the total fluxes. See text for parameter values. 
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Fig. 2.— High-energy spectrum of PKS 0208-512. The individual components are specified 
by the curve types as in Fig. 1. The parameters are ad = 0.5, © a = ©d = 0.5, e c = 10 ” 3 
7 i = 1, 72 = 10 6 , 7 ;, = 10 3 , pi = 2.1, p 2 = 3.1, T = 3 and 0 = 18°. 
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Fig. 3.— High-energy spectrum of PKS 0506-612. The individual components are specified 
by the curve types as in Fig. 1. The parameters are the same as in Figure 2 except V = 11 
and 0 = 9°. 



